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Abstract 

In this paper, we are going to discuss several approaches to solve the quadratic and linear 
simplicity constraints in the context of the canonical formulations of higher dimensional General 
Relativity and Supergravity developed in [U [2l |3] |U [5l |6j . Since the canonical quadratic simplicity 
constraint operators have been shown to be anomalous in any dimension D > 3 in [3J, non- 
standard methods have to be employed to avoid inconsistencies in the quantum theory. We 
show that one can choose a subset of quadratic simplicity constraint operators which are non- 
anomalous among themselves and allow for a natural unitary map of the spin networks in the 
kernel of these simplicity constraint operators to the SU(2)-based Ashtekar-Lewandowski Hilbert 
space in D = 3. The linear constraint operators on the other hand are non-anomalous by 
themselves, however their solution space will be shown to differ in D = 3 from the expected 
Ashtekar-Lewandowski Hilbert space. We comment on possible strategies to make a connection 
to the quadratic theory. Also, we comment on the relation of our proposals to existing work in the 
spin foam literature and how these works could be used in the canonical theory. We emphasise 
that many ideas developed in this paper are certainly incomplete and should be considered as 
suggestions for possible starting points for more satisfactory treatments in the future. 
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A Simple Irreps of SO(D + 1) and Square Integrable Functions on the Sphere 



1 Introduction 

In [H |2], gravity in any dimension D + 1 > 3 has been formulated as a gauge theory of S0(1, D) 
or of the compact group SO(D + 1), irrespective of the spacetime signature. The resulting theory 
has been obtained on two different routes, a Hamiltonian analysis of the Palatini action making 
use of the procedure of gauge unfixing^} and on the canonical side by an extension of the ADM 
phase space. The additional constraints appearing in this formulation, the simplicity constraints, 
are well known. They constrain bivectors to be simple, i.e. the antisymmetrised product of 
two vectors. Originally introduced in Plebanski's |10] formulation of General Relativity as a 
constrained BF theory in 3 + 1 dimensions, they have been generalised to arbitrary dimension in 
|llj and were considered in the context of Hamiltonian lattice gravity |12|I13]. Moreover, discrete 
versions of the simplicity constraints are a standard ingredient of the Spin Foam approaches to 
quantum gravity [21 [151 02], see JT7] [18] for reviews, and recently were also used in Group 
Field theory [191 E01 El] as well as on a simplicial phase space [221 123] , where also their algebra 
was calculated. Two different versions of simplicity constraints are considered in the literature, 
which are either quadratic or linear in the bivector fields. The quantum operators corresponding 
to the quadratic simplicity constraints have been found to be anomalous both in the covariant 
|24| as well as in the canonical picture |25l [3J. On the covariant side, this lead to one of the 

1 See [7| HIE] for original literature on gauge unfixing. 
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major points of critique about the Barrett-Crane model |14j : The anomalous constraints are 
imposed strongljj^J which may imply erroneous elimination of physical degrees of freedom |26j . 
This triggered the development of the new Spin Foam models [271 EE1 031 Ell EE EH] , i n which the 
quadratic simplicity constraints are replaced by linear simplicity constraints. The linear version 
of the constraints is slightly stronger than the quadratic constraints, since in 3 + 1 dimensions 
the topological solution is absent. The corresponding quantum operators are still anomalous 
(unless the Immirzi parameter takes the values 7 = iy^, where £ denotes the internal signature, 
or 7 = 00). Therefore, in the new models (parts of) the simplicity constraints are implemented 
weakly to account for the anomaly. Also, the newly developed U(iV) tools [3U1EI1E2] have been 
recently applied to solve the simplicity constraints [331 EH E2] . 

In this paper, we are first going to take an unbiased look at them from the canonical perspec- 
tive in the hope of finding new clues for how to implement the constraints correctly. Afterwards, 
we will compare our results to existing approaches from the Spin Foam literature and outline 
similarities and differences. We stress that will not arrive at the conclusion that a certain kind 
of imposition will be the correct one and thus further research, centered around consistency 
considerations and the classical limit, has to be performed to find a satisfactory treatment for 
the simplicity constraints. Of course, in the end an experiment will have to decide which imple- 
mentation, if any, will be the correct one. Since such experiments are missing up to now, the 
general guidelines are of course mathematical consistency of the approach, as well as compar- 
ison with the classical implementation of the simplicity constraints in D = 3, where the usual 
SU(2) Ashtekar variables exist. If a satisfactory implementation in D = 3 can be constructed, 
the hope would then be that this procedure has a natural generalisation to higher dimensions. 
Since parts of the very promising results developed from the Spin Foam literature are restricted 
to four dimensions, we will restrict ourselves to dimension independent treatments in the main 
part of this paper. 

The paper will be divided into three parts. We will begin with investigating the quadratic 
simplicity constraint operators which have been shown to be anomalous in [3]. It will be illus- 
trated that choosing a recoupling scheme for the intertwiner naturally leads to a maximal closing 
subset of simplicity constraint operators. Next, the solution to this subset will be shown to allow 
for a natural unitary map to the SU(2) based Ashtekar-Lewandowski Hilbert space in D = 3 
and we will finish the first part with several remarks on this quantisation procedure. In the 
section [3| we will analyse the strong implementation of the linear simplicity constraint operators 
since they are non-anomalous from start. The resulting intertwiner space will be shown to be 
one-dimensional, which is problematic because this forbids the construction of a natural map 
to the SU(2) based Ashtekar-Lewandowski Hilbert space. In contrast to the quadratic case, the 
linear simplicity constraint operators will be shown to be problematic when acting on edges. We 
will discuss several possibilities of how to resolve these problems and finally introduce a mixed 
quantisation, in which the linear simplicity constraints will be substituted by the quadratic con- 
straints plus a constraint ensuring the equality of the normals N 1 and n 1 ^). In section [ZJ we 
will compare our results to existing approaches from the Spin Foam literature. Finally, we will 
give a critical evaluation of our results and conclude in section [5] 

2 The Quadratic Simplicity Constraint Operators 

2.1 A Maximal Closing Subset of Vertex Constraints 

In our companion papers [U El El El El E]> a canonical connection formulation of (D + 1)- 
dimensional Lorentzian General Relativity was developed, using an SO(D + l)-connection A a u 

2 Strongly here means that the constraint operator annihilates physical states, C \ip) = V \ip) £ W P h V s 
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and its conjugate momentum ir aIJ as canonical variables. Here, a,b, . . . = 1,...,D are spatial 
tensorial indices and I, J, . . . = 0, . . . , D are Lie algebra indices in the fundamental representa- 
tion. A key input of the construction are the (quadratic) simplicity constraints 

which enforce, up to a topological sector present in D = 3, that ir aIJ ~ 2n^E a ^, where 
E aJ is an SO(D + 1) valued vector density, a so called hybrid vielbein, and n 1 is the unique 
(up to sign) normal defined by njE aI = 0. Fixing the time gauge n 1 = (1,0, ...,0), one 
arrives at the ADM (extended phase space) formulation of General Relativity with SO(D) gauge 
invariance, see [2] for details. The second class constraints which normally arise as stability 
conditions on the simplicity constraints are absent in our connection formulation, since they 
can be explicitly removed by the process of gauge unfixing after performing the Dirac analysis, 
see [2]. Essentially, they are gauge fixing conditions for the gauge transformations generated by 
the simplicity constraint, which change a certain part of the torsion of the A a jj. The square of 
this part of the torsion is included in a respective decomposition of the Palatini action and thus 
results in the second class partner for the simplicity constraint [2]. 

A quantisation of the simplicity constraint using loop quantum gravity methods results in 
a complicated operator, since TT aIJ becomes a flux operator which acts as the sum of all right 
invariant vector fields associated to the different edges at a vertex. In order to facilitate the 
treatment of this quantum constraint, it has been shown in [3] that the necessary and sufficient 
building blocks of the quadratic simplicity constraint operator acting on a vertex v are given by 

R[ijR e KL]f^ = Ve, e 'e{e"G E(j);v = b(e")}, (2.2) 

where Rjj is the right invariant vector field associated to the edge e, / 7 is the a cylindrical 
function defined on an adapted graph 7, e.g. a spin network, v is a vertex of 7, E{^) is the 
set of edges of 7 and 6(e) denotes the beginning of the edge e. The orientations of all edges 
are chosen such that they are outgoing of v. We note that these are exactly the off-diagonal 
simplicity constraints familiar from spin foam models, see e.g. |11| 124] . 

Since not all of these building blocks commute with each other, i.e. the ones sharing exactly 
one edge, we will have to resort to a non-standard procedure in order to avoid an anomaly in the 
quantum theory. The strong imposition of the above constraints, leading to the Barrett-Crane 
intertwiner [TJ] , was discussed in [11] . A master constraint formulation of the vertex simplicity 
constraint operator was proposed in [3] , however apart from providing a precise definition of the 
problem, this approach has not lead to a concrete solution up to now. 

In this paper, we are going to explore a different strategy for implementing the quadratic 
vertex simplicity constraint operators which is guided by two natural requirements: 

1. The imposition of the constraints should be non- anomalous. 

2. The imposition of the simplicity constraint operator in D = 3 should, at least on the 
kinematical level, lead to the same Hilbert space as the quantisation of the classical theory 
without a simplicity constraint. More precisely, there should exist a natural unitary map 
from the solution space of the quadratic simplicity constraint operators ^simple to the 
Ashtekar-Lewandowski Hilbert space %al in D = 3. 

The concept of gauge unfixing [7J [9] which was successfully used in order to derive the 
classical connection formulation of General Relativity [U [2] used in this paper was originally 
developed in the context of anomalous gauge theory, where it was observed that first class 
constraints can turn into second class constraints after quantisation [36, 37, 38, 391 HO]. This is 
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however precisely what is happening in our case: The classically Abelian simplicity constraints 
become a set of non-commuting operators due to the regularisation procedure used for the fluxes. 
The natural question arising is thus: How does a set of maximally commuting vertex simplicity 
constraint operators look like? 

Theorem 1. Given a N-valent vertex v G 7, the set 

€ IJKLM^eiRe 1 = • • • = e IJKLM^eifRe N = (2-3) 



€ IJKLM C^ei + R-ei + R-el 



e IJKLM [Ret + Re'2) (R^ + Re 2 L ) 

)(r? 1 l + r£ l + i£ l ) 



"e 3 ; 








6 IJKLM (-^ei 



I J 



+ ... + R 



IJ 



R« L + ... + R 



KL 

&N-2 







(2.4) 



generates a closed algebra of vertex simplicity constraint operators. Under the assumption that 
no linear combinations with different multi-indices M = M1M2 ■ ■ ■ Md—3 are allowed^ the set 
is maximal in the sense that adding new vertex constraint operators spoils closure. 

Proof. Closure can be checked by explicit calculation. In order to understand why the calculation 
works, recall that right invariant vector fields generate the Lie algebra so(D + 1) as |3j 



1 



Rh, RkL = o^,e' (VJKRlL + ViLRjK ~ VIkRjl ~ VJlR 



IK) 



(2.5) 



and thus infinitesimal rotations. The commutativity of (2.3) has been discussed in [3]. Further 



we see that every element of (2.4) operates on (2.3) as an infinitesimal rotation. The same is also 



true for the elements in (2.4): Taking the ordering from above, every constraint operates as an 



infinitesimal rotation on all constraints prior in the list. Since the commutator is antisymmetric 



in the exchange of its arguments, closure, i.e. commutativity up to constraints, of (2.4) follows 



To prove maximality of the set we will show that, having chosen a subset of simplicity 



constraints as given in (2.3) and (2.4), adding any other linear combination of the building 



blocks (2.2) spoils the closure of the algebra. To this end, we make the most general Ansatz 



E 



a 



y t IJKLM rL i j 



(2.6) 



l<i<j<N 



for an iV-valent vertex. Note that the diagonal terms (i = j) are proportional to (2.3) and 
therefore do not have to be taken into account in the above sum, and that Rjy 



can 



be dropped due to gauge invariance. Moreover, a« can be chosen such that for fixed j' not all 



(i < j') are equal. Otherwise, with on 



ij' 



a 



we find the term ctye IJKI j^R 



IJ 



R 



KL 



the sum, which can be expressed as a linear combination of (2.3) and (2.4) and therefore can be 
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A superposition of different multi-indices seems to be highly unnatural since an anomaly with the Gaufi 
constraint has to be expected. We are however currently not aware of a proof which excludes this possibility from 
the viewpoint of a maximal closing set. 
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dropped. Consider 



l e IJKLM R li R VZ '■> e abode {ai 3 Ri B R^ D + a 23 R 2 B R% D + ...)] 



nU f KL AB tjMN D 
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~ ^ ai 3 e IJKLM 1 
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m ^ - "2iJ £ljKLM K 2 J MnR\ ^ABCDE K j > 

i=3 



(2.7) 



of) simplicity building blocks (2.2 



where we dropped terms proportional to (2.3) in the first and in the second step. For a closing 
algebra, the right hand side of 

> 3) have to vanish separately 



■3- 



(In general, one could make use of gauge invariance to "mix" the contributions of different R 
However, in the case at hand this will produce terms containing Rn, which do not vanish if the 
contributions of different Rjs did not already vanish separately). 

We start with the case D = 3. The summands on the right hand sides of (2.7) are proportional 

to 



SfjW(R j )AB(R2y J (Ri)-c, 



\U i 



where we used the notation Sj"' 



■In 



n\ S\) 5 I, '..5 I 1 n , 

[J l J2 Jn\ 



To show that this expression can not be 



rewritten as a linear combination of the of building blocks (2.2) we antisymmetrise the indices 



[ABU], [ABKC] and [IJKC] and find in each case that the result is zero. 
For D > 3, the summands are proportional to 



zABC_E 



(R j ) AB (R 2 ) IJ (R 1 ) K c 



(2.9) 



Whatever multi-index E we might have chosen in the Ansatz (2.6), we can always restrict 



attention to those simplicity constraints in the maximal set which have the same multi-index 
M = E. Then, the same calculation as in the case of D = 3 shows that the antisymmetrisations 
of the indices [ABU], [ABKC] and [IJKC] vanish. 

Therefore, the only possibilities are (a) the trivial solution ay = a 2 j = or (b) a±j = 



a 2j(7^ 0), which implies that the terms on the right hand side of (2.7) are a rotated version of 



1 j 



r>KL 
K 2 ■ 



e IJKLM R 

have ai3 = 023 
we can define aL 



The second option (b) is, for j = 3, 



excluded by our choice of 



and we must 

0. Next, consider j = 4 and suppose we have 014 = 024 := a' 7^ 0. Then, 
34 := 034 — a' and find the terms o! e IJKJ ^R^ 23) R^ L + 0£' M e IJKL jjRl J Rf L in 
(2.6). The first term again is already in the chosen set, which implies we can set 0:14 = a 2 4 = 
w.l.o.g. by changing 034 — > a 34 (We will drop the prime in the following 
generalises to j > 4, and we have w.l.o.g. ay = a 2 j = (3 < j < N). 

Suppose we have calculated the commutators of £jjklm R i 
and found that for closure, we need ctjj 



This immediately 



(i = 2,.. 

for 1 < i < n and i < j < N . Then, 



?IJ ,Ri L 



,n) with (2.6) 



N-l 
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(2.10) 
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which, by the reasoning above, again is not a linear combination of any simplicity building blocks 
for any choice of ar n+ uj, and therefore only the trivial solution ai n -\-\)j = (n + 1 < j < N) 
leads to closure of the algebra. □ 

2.2 The Solution Space of the Maximal Closing Subset 



In order to interpret this set of constraints recall from [3] that the constraints in (2.3) are the 
same as the diagonal simplicity constraints acting on edges of 7 and can be solved by demanding 



the edge representations to be simple. The remaining constraints (2.4) can be interpreted as 



specifying a recoupling scheme for the intertwiner 1 at v: Couple the representations on e\ and 
e2, then couple this representation to e^, and so forth, see fig. [TJ We call the intermediate virtual 
edges ei2, ei23, • • • and denote the highest weights of the representations thereon by A12, A123, . . . 
Since we can use gauge invariance at all the intermediate intertwiners in the recoupling scheme, 
e.g., R ei + R e2 = R ei2 , we have 



_ (n 1J _i_ j?iJ\ ivi KL _i_ u KL \ — a —P IJ p 

IJKLM \ ri e 1 T -n-e 2 J \ rL e 1 ' ^2 J ~ t IJKLM ri e 1 2 lx , 

and thus that the representation on e\2 has to be simple, i.e. 



l ei2 



A 



12 



(Ai2,0,...,0) A12 = 0,1,2,.. 



(2.11) 



(2.12) 



Using the same procedure, all intermediate representations are required to be simple and the 
intertwiner is labeled by N — 3 "spins" Aj G No- We call an intertwiner where all internal lines 
are labeled with simple representations simple. 




Figure 1: Recoupling scheme corresponding to the subset of quadratic vertex simplicity con- 



straint operators (2.4). 



Denote byl^^ the set of iV-valent SU(2) intertwiners and by j;Sp™( d + 1 ) ^ e se ^ f simple iV- 
valent Spin(D + 1) intertwiners. Recalling that an iV-valent SU(2) intertwiner can be expressed 
in the same recoupling basis and calling the intermediate spins ji, we see that the map 

p . T Spin(D+l) T SU(2) 
b ■ 1 s,N N 

\h ^ ji (2.13) 

is unitary (with respect the scalar products induced by the respective Ashtekar-Lewandowski 
measures, see [3]). The motivation for the factor 1/2 comes from the fact that A = (1,0) in 
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D = 3 corresponds to the familiar j + = j- = 1/2 and the area spacings of the SO(4) and the 
SU(2) based theories agree using this identification, cf [3]. 



2.3 Remarks 

1. Since the choice of the maximal closing subset of the simplicity constraint operators is 
arbitrary, no recoupling basis is preferred a priori. On the SU(2) level, a change in the 
recoupling scheme amounts to a change of basis in the intertwiner space and therefore poses 
no problems. On the level of simple Spin(D + 1) representations however, a choice in the 
recoupling scheme affects the property "simple" , since the non-commutativity of constraint 
operators belonging to different recoupling schemes means that kinematical states cannot 
have the property simple in both schemes. 

2. There exist recoupling schemes which are not included in the above procedure, e.g., take 
N = 6 and the constraints eR^Ru = £-^34-^34 = eR^R^ = and couple the three re- 
sulting simple representations. The theorem should however generalise to those additional 
recoupling schemes. 

3. It is doubtful if the action of the Hamiltonian constraint leaves the space of simple inter- 
twiners in a certain recoupling scheme invariant. To avoid this problem, one could use a 
projector on the space of simple intertwiners in a certain recoupling scheme to restrict the 
Hamiltonian constraint on this subspace and average later on over the different recoupling 
schemes if they turn out to yield different results. The possible drawbacks of such a pro- 
cedure are however presently unclear to the authors and we refer to further research. The 
construction of such a projector can be seen as a quantum analogue of the gauge unfixing 
process familiar from our companion paper [2]. A possible strategy to find a Hamiltonian 
constraint operator which leaves the solution space of a first class subset invariant is to 
construct a gauge unfixing projector which adds vertex simplicity constraints which are 
not in the first class subset to the Hamiltonian constraint such that it commutes with the 
first class subset. 

4. It would be interesting to check whether the dropped constraints are automatically solved 
in the weak operator topology (matrix elements with respect to solutions to the maximal 



5. The imposition of the constraints can be stated as the search for the joint kernel of a 
maximal set of commuting generalised area operators 



Notice, however, that for D > 3 these generalised area operators, just as the simplicity 
constraints, are not gauge invariant while in D = 3 they are 

6. In D = 3 we have the following special situation: 

We have two classically equivalent extensions of the ADM phase at our disposal whose 
respective symplectic reduction reproduces the ADM phase space. One of them is the 
Ashtekar-Barbero-Immirzi connection formulation in terms of the gauge group SU(2) with 
additional SU(2) Gaufi constraint next to spatial diffeomorphism and Hamiltonian con- 
straint, and the other is our connection formulation in terms of SO (4) with additional 
SO(4) Gaufi constraint and simplicity constraint. Both formulations are classically com- 
pletely equivalent and thus one should expect that also the quantum theories are equivalent 



subset). 




(2.14) 
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in the sense that they have the same semiclassical limit. Let us ask a stronger condition, 
namely that the joint kernel of SO(4) Gaufi and simplicity constraint of the SO(4) theory 
is unitarily equivalent to the kernel of the SU(2) Gaufi constraint of the SU(2) theory. To 
investigate this first from the classical perspective, we split the SO (4) connection and its 
conjugate momentum (A IJ ,mj) into self-dual and anti-selfdual parts A^,^) which then 
turn out to be conjugate pairs again. It is easy to see that the SO(4) Gaufi constraint Gu 
splits into two SU(2) Gaufi constraints , one involving only self-dual variables and the 
other only anti-selfdual ones which therefore mutually commute as one would expect. The 
SO(4) Gaufi constraint now asks for separate SU(2) gauge invariance for these two sec- 
tors. Thus a quantisation in the Ashtekar-Isham-Lewandowski representation would yield 
a kinematical Hilbert space with an orthonormal basis T^<g>T~_ where S± are usual SU(2) 
invariant spin networks. The simplicity constraint, which in D = 3 is Gaufi invariant and 
can be imposed after solving the Gaufi constraint, from classical perspective asks that the 
double density inverse metrics q± = 7r?±7<"jy-£ are identical. This is classically equivalent 
to the statement that corresponding area functions Ar±(S) are identical for every S. The 
corresponding statement in the quantum theory is, however, again anomalous because it 
is well known that area operators do not commute with each other. On the other hand, 
neglecting this complication for a moment, it is clear that the quantum constraint can only 
be satisfied on vectors of the form T+_ (g) T~_ for all S if s+, s_ share the same graph and 
SU(2) representations on the edges because if S cuts a single edge transversally then the 
area operator is diagonal with an eigenvalue oc + 1) and we can always arrange such 

an intersection situation by choosing suitable S. By a similar argument one can show that 
the intertwiners at the edges have to be the same. But this is only a sufficient condition 
because in a sense there are too many quantum simplicity constraints due to the anomaly. 
However, the discussion suggests that the joint kernel of both SO(4) and simplicity con- 
straint is the closed linear span of vectors of the form T+ ® T~ for the same spin network 
s = s+ = S- . The desired unitary map between the Hilbert spaces would therefore simply 
be T s H- T+ ® T~ . 

This can be justified abstractly as follows: From all possible area operators pick a maximal 
commuting subset Ar^ using the axiom of choice (i.e. pick a corresponding maximal set 
of surfaces S a ). We may construct an adapted orthonormal basis T^ 1 diagonalising all of 
thenj^Jsuch that Ar^Tj^ = X a T^. Now the constraint 

Ar+ <g> 1 = 1 ® Ar~ 

can be solved on vectors by demanding A + = A_ . The desired unitary map 
would then be T\ i— y (g) . Thus the question boils down to asking whether a maximal 
closing subset can be chosen such that the eigenvalues A are just the spin networks s. We 
leave this to future research. 

7. the afore mentioned split into selfdual and anti-selfdual sector is meaningless and 

we must stick with the dimension independent scheme outlined above. An astonishing fea- 
ture of this scheme is that after the proposed implementation of the simplicity constraints, 
the size of the kinematical Hilbert space is the same for all dimensions D > 3/ By "size", 
we mean that the spin networks are labelled by the same sets of quantum numbers on the 
graphs. Of course, before imposing the spatial diffeomorphism constraint these graphs are 

4 If the maximal set still separates the points of the classical configurations space, this should leave no room 
for degeneracies, that is the \ a completely specify the eigenvector. We will assume this to be the case for the 
following argument. 
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embedded into spatial slices of different dimension and thus provide different amounts of 
degrees of freedom. However, after implementation of the diffeomorphism constraint, most 
of the embedding information will be lost and the graphs can be treated almost as abstract 
combinatorial objects. Let us neglect here, for the sake of the argument, the possibility 
of certain remaining moduli, depending on the amount of diffeomorphism invariance that 
one imposes, which could a priori be different in different dimensions. In the case that the 
proposed quantisation would turn out to be correct, that is, allow for the correct semiclas- 
sical limit, this would mean that the dimensionality of space would be an emergent concept 
dictated by the choice of semiclassical states which provide the necessary embedding infor- 
mation. A possible caveat to this argument is the remaining Hamiltonian constraint and 
the algebra of Dirac observables which critically depend on the dimension (for instance 
through the volume operator or dimension dependent coefficients, see jU |2]) and which 
could require to delete different amounts of degrees of freedom depending on the dimension. 

This idea of dimension emergence is not new in the field of quantum gravity, however, 
it is interesting to possibly see here a concrete technical realisation which appears to be 
forced on us by demanding anomaly freedom of the simplicity constraint operators. Of 
course, these speculations should be taken with great care: The number of degrees of free- 
dom of the classical theory certainly does strongly depend on the dimension and therefore 
the speculation of dimension emergence could fail exactly when we try to construct the 
semiclassical sector with the solutions to the simplicity constraints advertised above. This 
would mean that our scheme is wrong. On the other hand, there are indications [41J that 
the semiclassical sector of the LQG Hilbert space already in D = 3 is entirely described in 
terms of 6-valent vertices. Therefore, the higher valent graphs which in D = 3 could cor- 
respond to pure quantum degrees of freedom, could account for the semiclassical degrees 
of freedom of higher dimensional General Relativity. Since there is no upper limit to the 
valence of a graph, this would mean that already the D = 3 theory contains all higher 
dimensional theories! 

Obviously, this puzzle asks for thorough investigation in future research. 

8. The discussion reveals that we should compare the amount of degrees of freedom that 
the classical and the quantum simplicity constraint removes. This is a difficult subject, 
because there is no well defined scheme that attributes quantum to classical degrees of 
freedom unless the Hilbert space takes the form of a tensor product, where each factor 
corresponds to precisely one of the classical configuration degrees of freedom. The follow- 
ing "counting" therefore is highly heuristic and speculative: 

In the case D = 3, the classical simplicity constraints remove 6 degrees of freedom from 
the constraint surface per point on the spatial slice. In order to count the quantum degrees 
of freedom that are removed by the quantum simplicity constraint when acting on a spin 
network function, we make the following, admittedly naive analogy: 

We attribute to a point on the spatial slice an N- valent vertex v of the underlying graph 
7 which is attributed to the spatial slice. This point is equipped with degrees of freedom 
labelled by edge representations and the intertwiner. Every edge incident at v is shared 
by exactly one other vertex (or returns to v which however does not change the result). 
Therefore, only half of the degrees of freedom of an edge can be attributed to one vertex. 
We take as edge degrees of freedom the L^^J Casimir eigenvalues of SO(D + 1) labelling 
the irreducible representation. The edge simplicity constraint removes all but one of these 
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Casimir eigenvalues, thus per edge [— 5— J edge degrees of freedom are removed. Further, a 
gauge invariant intertwiner is labelled by a recoupling scheme involving N — 3 irreducible 
representations not fixed by the irreducible representations carried by the edges adjacent 
to the vertex in question, which are fully attributed to the vertex (there are N — 2 virtual 
edges coming from coupling 1,2 then 3 etc. until N but the last one is fixed due to gauge 
invariance). We take as vertex degrees of freedom these N — 3 irreducible representa- 
tions each of which is labelled again by [^x^J Casimir eigenvalues. The vertex simplicity 
constraint again deletes all but one of these eigenvalues, thus it removes (N — 
quantum degrees of freedom. We conclude that the quantum simplicity constraint removes 

(A r -3 + f)L^ 1 J (2-15) 

quantum degrees of freedom per point (iV-valent vertex) where N—3 accounts for the vertex 
and N/2 for the TV" edges counted with half weight as argued above. This is to be compared 
with the classical simplicity constraint which removes D 2 (D — l)/2 — D degrees of freedom 
per point. Requiring equality we see that vertices of a definitive valence Nd are preferred 
in D spatial dimensions which for large D grows quadratically with D. Specifically for 
D = 3 we find N3 = 6. Thus, our naive counting astonishingly yields the same preference 
for 6-valent graphs in D = 3 as has been obtained in j^T] by completely different methods. 
From the analysis of |41j . it transpires that N3 = 6 has an entirely geometric origin and 
one thus would rather expect Np = 2D (hypercubulations) and this may indicate that our 
counting is incorrect. 



3 The Linear Simplicity Constraint Operators 
3.1 Regularisation and Anomaly Freedom 

In [5] , the connection formulation sketched at the beginning of the previous chapter was altered in 
that it contains linear simplicity constraints e ijKLM^ l7T<lJK ~ ^ an< ^ an independent normal N 1 
as phase space variables. The normal Poisson-commutes with both the connection A a jj and its 
momentum 7r aIJ and has its own canonical momentum Pj. The necessity for this independent 
normal did not stem from the anomaly encountered when looking at the quadratic quantum 
simplicity constraints, but from the observation that it was needed to extend the connection 
formulation to higher dimensional supergravities. 

Since the linear simplicity constraint is a vector of density weight one, it is most naturally 
smeared over [D — l)-dimensional surfaces. The regularisation of the objects 

S b (S):= f b Ljl (x)e IJKm N I (x)7r aJK (x)e abl ... bD ^ 1 dx^ A...Adx b °~\ (3.1) 
Js 

where S b denotes the linear simplicity constraint, S a D — 1-surface, and b LM an arbitrary 
semianalytic smearing function of compact support, therefore is completely analogous to the 
case of flux vector fields. The corresponding quantum operator 

§\S)f = f**(S)f = P ; s Y*f(S)f ys = P ; s £ e (e, S)e IJKlM b L ' M (b(e))N I (b{e))R J e K f ls 

(3.2) 

has to annihilate physical states for all surfaces S C a and all semianalytic functions b IM of 
compact support, where p 7 denotes the cylindrical projection and 75 is a graph adapted to the 
surface S. Since we can always choose surfaces which intersect a given graph only in one point, 
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this implies that the constraint has to vanish when acting on single points of a given graph. In 
[3], it has been shown that the right invariant vector fields actually are in the linear span of the 
flux vector fields. Therefore, it is necessary and sufficient to demand that 



e IJKLM b 



LM 



(b(e)) fi'{b(e)) Ri K ■ f 7 = 



(3.3) 



for all points of 7 (which can be be seen as the beginning point of edges by suitably subdividing 
and inverting edges). Since N 1 acts by multiplication and commutes with the right invariant 
vector fields, see [5] for details, the condition is equivalent tc[^] 



R IJ ■ / 7 = 0, 



(3.4) 



i.e. the generators of rotations stabilising N 1 have to annihilate physical states. Before imposing 
these conditions on the quantum states, we have to consider the possibility of an anomaly. 
Classically and before using the singular smearing of holonomies and fluxes, both, the linear 
and the quadratic simplicity constraints are Poisson self-commuting. The quadratic constraint 
is known to be anomalous both in the Spin Foam [24J as well as in the canonical picture \25\ [3] 
and thus should not be imposed strongly. Also the linear simplicity constraint is anomalous 
when using a non-zero Immirzi parameter (at least if 7 7^ 1 in the Euclidean theory. But 
7 = 1 is ill-defined for SO (4), see e.g., 02] )• Surprisingly, in the case at hand and without an 
Immirzi parameter in four dimensions, we do not find an anomaly. However that is just because 
the generators of rotations stabilising N 1 form a closed subalgebra! Direct calculation yields, 
choosing (without loss of generality) 755' to be a graph adapted to both surfaces S, S', 



Sj ss ,( S )iSj S s'( S ") Ass' 



...K e , ^ • • • U e' 

e£7ss' e '£7ss' 



Ass' ~~ E ' " [ R i J i R e B ] f-TsS' 
ee 7ss' 



E*I «J -B f KL CD 

■ ■■V kV lV cV d j 



AIN 



nMN f 



e S7ss' 



= £ . . . rfW^crfn {rj L ^8 D \ M 8 N f) R?» / 7ss , 

e S7sS' 
e S7ss' 

where the operator in the last line is in the linear span of the vector fields S b (S). The classical 
constraint algebra is not reproduced exactly (the commutator does not vanish identically), but 
the algebra of quantum simplicity constraints closes, they are of the first class. Therefore, strong 
imposition of the quantum constraints does make mathematical sense. 

Note that up to now, we did not solve the Gaufl constraint. The quantum constraint algebra 
of the simplicity and the GauB constraint can easily be calculated and reproduces the classical 



5 Use the decomposi tion of Xu into its rotational (Xu := r)f rjj Xkl) and "boost" parts (Xi := —C,N J Xu) 

(N'N J and C = 1 for SO(D + 1) and C = -1 for SO(l,D) as 



with respect to N 1 in (3.3 1, where r/ IJ = r) IJ 



internal gauge groups. It follows that 77 Nj — 0. 
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result 



s b (s),G AB [A AB ] P ;j JS 



? AB 

Iff 



IS 



E £ ( e '' S^^b^N 1 ^/, A AB (v) £ + iJ 

e'6E(7s),'y=&(e') \ee^(7 S ),t)=6(e) 
eEE(js),v=b(e) 

^A^Ct;)^^^) E e(e, 5) (^(«)2^^ B + ^ A N B \v)R J e K ) / 7s 

ee£(7 S ),-u=&(e) 

=^- A - fo )(5) P ; s / 7s . (3.6) 

where we used := \ (j^ A gfj^ ~ ^al^) • ^ follows that the simplicity constraint operator 
does not preserve the Gaufi invariant subspace (in other words, as in the classical theory, the 
Gaufi constraint does not generate an ideal in the constraint algebra). This implies that the joint 
kernel of both Gaufi and simplicity constraint must be a proper subspace of the Gaufi invariant 
subspace. It is therefore most convenient to look for the joint kernel on in the kinematical (non 
Gaufi invariant) Hilbert space. 

3.2 Solution on the Vertices 

Consider a slight modification of the usual gauge- variant spin network functions, where the 
intertwiners i v = i v (N) are square integrable functions of N 1 . Let v be a vertex of 7 and 
ei, . . . , e n the edges of 7 incident at v, where all orientations are chosen such that the edges are 
all outgoing at v. Then we can write the modified spin network functions 



i=i 1 1 

= tx(i v (N)-®? =1 n ei (h ei (A))-M v ), (3.7) 

where M v contracts the indices corresponding to the endpoints of the edges e, and represents the 
rest of the graph 7. These states span the combined Hilbert space for the normal field and the 
connection Ht = %grav ®T~Ln (cf. [5]) and they will prove convenient for solving the simplicity 
constraints. Choose the surface S' such that it intersects a given graph 7' only in the v ertex 



v' G 7'. The action of S b (S') on the vertex v' of a spin network T , ^(A,N) implies with (3.4) 
that 

S b (S%,T YX j(A,N) = 
tr((i„(iV)r^) -®? =1 n Bi {h e M))-M v ) = Ve at v', (3.8) 

where here denote the generators of SO(D + 1) in the representation 7r; e of the edge e 
and the bar again denotes the restriction to rotational components (w.r.t. N 1 ). The above 
equation implies that the intertwiner i v , seen as a vector transforming in the representation 
7f/ e dual to 7i7 e of the edge e, has to be invariant under the SO(D)n subgroup which stabilises 
the N 1 . By definition [33], the only representations of SO(D + 1) which have in their space 
nonzero vectors which are invariant under a SO(-D) subgroup are of the representations of class 
one (cf. also appendix Kl), and they exactly coincide with the simple representations used in 
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Spin Foams It is easy to see that the dual representations (in the sense of group theory) 
of simple representations are simple representations. Therefore, all edges must be labelled by 
simple representations of SO(D + 1). Moreover, SO(D) is a massive subgroup of SO(D + 1) 
|43j . so that the (unit) invariant vector £i e (N) in the representation 7fj e is unique, which implies 
that the allowed intertwiners i v (N) are given by the tensor product of the invariant vectors of 
all n edges and potentially an additional square integrable function F V (N), i v (N) = £/ e (N) ® 
... (8) £i e (N) ® F V (N). Going over to normalised gauge invariant spin network functions implies 
that F V (N) = 1, and the resulting intertwiner space solving the simplicity and GauB constraint 
becomes one-dimensional, spanned by I V (N) := £j e (TV) ® ... (g) £i e (N). We will call these 
intertwiners and vertices coloured by them linear-simple. For an instructive example of the 
linear-simple intertwiners, consider the defining representation (which is simple since the highest 
weight vector is A = (1,0, ...,0), cf. appendix [A]) . The unit vector invariant under rotations 
(w.r.t. N 1 ) is given by N 1 and for edges in the defining representation incoming at v we simply 
contract h e J Nj. If the constraint is acting on an interior point of an analytic edge, this point 
can be considered as a trivial two-valent vertex and the above result applies. Since this has 
to be true for all surfaces, a spin network function solving the constraint would need to have 
linear-simple intertwiners at every point of its graph 7, i.e. at infinitely many points, which is in 
conflict with the definition of cylindrical functions (cf. [44] ) . In the next section, we comment 
on a possibility of how to implement this idea. 

3.3 Edge Constraints 

As noted above, the imposition of the linear simplicity constraint operators acting on edges is 
problematic, because it does not, as one might have expected, single out simple representations, 
but demand that at every point where it acts, there should be a linear-simple intertwiner. The 
problem with this type of solution is that all intertwiners, even trivial intertwiners at all interior 
points of edges, have to be linear-simple, which is however in conflict with the definition of a 
cylindrical function, in other words, there would be no holonomies left in a spin network because 
every point would be a A-dependent vertex. 

It could be possible to resolve this issue using a rigging map construction [35j HU H7] of the 
type 

r)(T rr A [T , f, p *] := lim C (pu, T 7 , Ty) (t p Z t T , r, r, Tl ) , (3.9) 

where *P 7 is the set of finite point sets p of a graph 7, p = {{xi}f =l \xi G jyi,N < 00}. ^3 7 
is partially ordered by inclusion, q >z p if p is a subset of q, so that the limit is meant in the 
sense of net convergence with respect to *B 7 . By the prescription T P7 .^ - we mean the projection 

of j-j-^ r onto linear-simple intertwiners at every point in p and C (p 7 , Ty, Ty) is a numerical 
factor. Assuming this to work, consider any surface S intersecting 7'. We (heuristically) find 

^[S^Ty) = lim C (p 7 , T 7 , Ty) (l*, 5 6 (S)T y ) 

= lim C(p„T„Ty)([S b (S^Ti;,T Y ) 

qj 7 9p 7 — ^00 \ 'I kin 

lim C(p 7 ,T 7 ,Ty) (s b (S)TP,T 1 ,) =0, (3.10) 

s P 7 9p 7 -s>oo x ' ' ' ' \ ' ' / kin 

since the intersection points of S with 7 will eventually be in p 7 and S b (S) is self-adjoint. 

We were however not able to find such a rigging map with satisfactory properties. It is 
especially difficult to handle observables with respect to the linear simplicity constraint and to 
implement the requirement, that the rigging map has to commute with observables. It therefore 
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seems plausible to look for non-standard quantisation schemes for the linear simplicity constraint 
operators, at least when acting on edges. Comparison with the quadratic simplicity constraint 
suggests that also the linear constraint should enforce simple representations on the edges, see 
the following remarks as well as section 3.5 for ideas on how to reach this goal. 



3.4 Remarks 

The intertwiner space at each vertex is one-dimensional and thus the strong solution of the un- 
altered linear simplicity constraint operator contrasts the quantisation of the classically imposed 
simplicity constraint at first sight. A few remarks are appropriate: 

1. One could argue that the intertwiner space at a vertex v is infinite-dimensional by taking 
into account holonomies along edges e' originating at v and ending in a 1-valent vertex 
v' . Since e and v' are assigned in a unique fashion to v if the valence of v is at least 2, 
we can consider the set {v,e',v'} as a new "non-local" intertwiner. Since we can label e' 
with an arbitrary simple representation, we get an infinite set of intertwiners which are 
orthogonal in the above scalar product. This interpretation however does not mimic the 
classical imposition of the simplicity constraints or the above imposition of the quadratic 
simplicity constraint operators. 

2. The main difference between the formulation of the theory with quadratic and linear 
simplicity constraint respectively is the appearance of the additional normal field sector 
in the linear case. Thus one could expect that one would recover the quadratic simplicity 
constraint formulation by ad hoc averaging the solutions of the linear constraint over the 
normal field dependence with the probability measure vm defined in [5]. Indeed, if one 
does so, then one recovers the solutions to the quadratic simplicity constraints in terms 
of the Barrett-Crane intertwiners in D = 3 and higher dimensional analogs thereof as has 
been shown long ago by Freidel, Krasnov, and Puzio [llj. A similar observation has been 
made in |48| . Such an average also deletes the solutions with "open ends" of the previous 
remark by an appeal to Schur's lemma. Since after such an average the N dependence 
of all solutions disappears, we can drop the integral in the kinematical inner product 
since /xtv is a probability measure. The resulting effective physical scalar product would 
then be the Ashtekar-Lewandowski scalar product of the theory between the solutions to 
the quadratic simplicity constraints. Such an averaging would also help with the solution 
of the edge constraints, since a 2-valent linear-simple intertwiner is averaged as 

/ dv{N)lf{N)^{N) = ^-5^, (3.11) 
Js D d n 

thus yielding a projector on simple representations. 

3. It can be easily checked that the volume operator as defined in [3], and therefore also more 
general operators like the Hamiltonian constraint, do not leave the solution space to the 
linear (vertex) simplicity constraints invariant. A possible cure would be to introduce a 
projector Vs on the solution space and redefine the volume operator as V := VsVVs- Such 
procedures are however questionable on the general ground that anomalies can always be 
removed by projectors. 

4. If one accepts the usage of the projector Vs, calculations involving the volume operator 
simplify tremendously since the intertwiner space is one-dimensional. We will give a few 
examples which can be calculated by hand in a few lines, restricting ourselves to the 
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defining representation of SO(D + 1), where the SO(-D)jv invariant unit vector is given by 
N 1 . 

Having direct access to N 1 , one can base the quantisation of the volume operator on the 
classical expression 



det q 



D\ 



.a D 



l 

D-l 



(3.12) 



In the case D + l uneven, this choice is much easier than the expression quantised in [3]. In 
the case -D + l even, the above choice is of the same complexit jP] as the one in [3], but leads 
to a formula applicable in any dimension and therefore, for us, is favoured. Proceeding as 
in [3], we obtain for the volume operator 



V{R) 
V(p) 



d u p\de%{q)(p)\ 



V, 



v ,7 



Qe 1 ,...,e D 



R 

D 



D 



d n P V(p) 1 , 



7' 



I 



^2 s(e 1 ,...,e D )q eq: ... :eD 

ei,...,e D eE(-y), e\<T\...C\e D =v 



(D-l) 



(3.13) 
(3.14) 

(3.15) 
(3.16) 



Note that the operator q ei ... eD is built from D right invariant vector fields. Since these 



are antisymmetric, q^ x er} 



In the case at hand, we have to use the 



projectors Vs to project on the allowed one-dimensional intertwiner space, the operator 
VsqPs therefore has to vanish for the case D + l even (an antisymmetric matrix on a one- 
dimensional space is equal to 0). However, the volume operator depends on q 2 , and Vsq^Vs 
actually is a non-zero operator in any dimension, though trivially diagonal. Therefore, also 
V is diagonal. 

The simplest non-trivial calculation involves a D-valent non-degenerate (i.e. no three 
tangents to edges at v lie in the same plane) vertex v where all edges are labelled by the 
defining representation of SO(D + 1) and thus the unique intertwiner which we will denote 
by |iV Al . . .N Ad ). We find 



q ei ,...,e D \N Al ...N A ») 



1 



S ( ei ,...,e D )(--) \Nje IAl '" Al ' 



D 



q eu ..., eD \Nje IA - A -) = s( ei ,...,e D )(iy D\\N A K..N A °), 



V v \N Al . ..N Ad ) 



D \ 2(D-1) 

D\ 



\N Al ...N Ad ) 



(3.17) 



i.e. for those special vertices, the volume operator preserves the simple vertices. For 
vertices of higher valence and/or other representations, we need to use the projectors. Of 



6 Up to (N) D+1 , but in the chosen representation N acts by multiplication and therefore is less problematic 
than additional powers of right invariant vector fields. 
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special interest are the vertices of valence D + l (triangulation) and 2D, where every edge 
has exactly one partner which is its analytic continuation through v (cubulation). We find 



V v \N Al ...iV AD + 1 } = - (D + 1)M \N Al ...N^ 1 ) , 




2(D-1) 



V v \N Al ...N A2D , cubic) = - (D)! \N Al ... N A ' 2D , cubic) . (3.1 



The dimensionality of the spatial slice now appears as a quantum number like the spins 
labelling the representations on the edges and it could be interesting to consider a large 
dimension limit in the spirit of the large N limit in QCD. 

5. When introducing an Immirzi parameter in D = 3 |2j, i.e. using the linear constraint 
ei.jKLN J TT aKL « while having {A aU {x), ^TT bKL {y)} = 2S b a 5^ 5^ D \x-y) with M7r a/J = 

TT aIJ + \/{2 1 )e IJKL K a KL, the linear simplicity constraint operators become anomalous 
unless 7 = +\/C> the (anti)self-dual case, which ho weve r results in non-invertibility of 
the prescription W. Repeating the steps in section 



3.1 



we find that these anomalous 

KL 



constraints require e IJKL N I (R^ L - 1/ (2 1 )e KLMN R^ IN ) • / 7 = 0. Since cuklN'^R, 
I / (2j)e K LM N R^ 1 N ) do not generate a subgroup, the constraint can not be satisfied strongly 
if the edge e transforms in an irreducible representation of SO(D + 1) (by definition, the 
representation space does not contain an invariant vector). 



In order to figure out the "correct" quantisation, one can try, in analogy to the strategy for 
the quadratic simplicity constraints, to weaken the imposition of the constraints at the quantum 
level. The basic difference between the linear and the quadratic simplicity constraints is that 
the time normal N is left arbitrary in the quadratic case and fixed in the linear case. In order 
to loose this dependence in the linear case, one could average over all N 1 at each point in a, 
which however leads to the Barrett-Crane intertwiners as described above. In analogy to the 
quadratic constraints, we could choose the subset 

e IJKLM NJ ( R e! L + R e 2 ) = 
e IJKLM NJ ( R fx L + R e 2 L + R e 3 ) = 

CijklmN J (< L + • • . + Rf N L _ 2 ) = (3.19) 

for each iV-valent vertex plus the edge constraints. As above, the choice of the subset specifies 
a recoupling scheme and the imposition of the constraints leads to the contraction of the virtual 
edges and virtual intertwiners of the recoupling scheme with the SO(D)jy-invariant vectors 
£l e .(N) and their complex conjugates ^ e .(-ZV), see fig. pi Gauge invariance can still be used 
at each (virtual) vertex in this calculation in the form j Rei = 0, which is sufficient since 
only R ei appears in the linear simplicity constraints. If we now integrate over each pair of 
£i e _ (N) "generated" by the elements of the proposed subset of the simplicity constraint operators 
separately, we obtain projectors on simple representations for each of the virtual edges in the 
recoupling scheme. The integration over N 1 for the edge constraints yields projectors on simple 
representations in the same manner. Finally, we obtain the simple intertwiners of the quadratic 
operators in addition to solutions where incoming edges are contracted with SO(D)Ar-invariant 
vectors £i e .(N). A few remarks are appropriate: 
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Figure 2: Recoupling scheme corresponding to the subset of linear vertex simplicity constraint 



operators (3.19). 



6. Although this procedure yields a promising result, it contains several non-standard and 
ah-hoc steps which have to be justified. One could argue that the "correct" quantisation 
of the linear and quadratic simplicity constraints should give the same quantum theory, 
however, as is well known, classically equivalent theories result in general in non-equivalent 
quantum theories, which nevertheless can have the same classical limit. 

7. It is unclear how to proceed with "integrating out" N 1 in the general case. For the vacuum 
theory, integration over every point in a gives the Barrett-Crane intertwiner for the edges 
contracted with SO(-D)jv-invariant vectors. This type of integration would also get rid of 
the 1-valent vertices and thus allow for a natural unitary map to the quadratic solutions 
as already mentioned above. 

8. When introducing fermions, there is the possibility for non-trivial gauge-invariant functions 
of N 1 at the vertices which immediately results in the question of how to integrate out 
this iV^-dependence. Next to including those N 1 in the above integration or to integrate 
out the remaining N 1 separately, one could transfer this integration back into the scalar 
product. Since the authors are presently not aware of an obvious way to decide about 
these issues, we will leave them for further research. 



3.5 Mixed Quantisation 

Since the implementation of the quadratic simplicity constraints described above yields a more 
promising result than the implementation of the linear constraints, we can try to perform a 
mixed quantisation by noting that we can classically express the linear constraints for even D 
in the form 

\ziJKLMK aIJ K bKL « 0, iV'-n'^O. (3.20) 

The phase space extension derived in [5] remains valid when interchanging the linear simplicity 
constraint for the above constraints. The reason for restricting D to be even is that we have 
an explicit expression for n J (-7r), see [H [2]. Since a quantisation of n 1 ^) will most likely not 
commute with the Hamiltonian constraint operator, we resort to a master constraint. Note that 
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the expression 

M , _ ((N 1 ~ »fo)Vg^ Sjj ((jV J - ^(tt))^- 1 ) 
which is the densitised square of N 1 — n / (-7r), can be quantised as 



M' N = 2^f- 2D (y/\V I V I \ - NjV 1 ), (3.22) 

when using a suitable factor ordering, where a quantisation of y/q 3 2D is described in [3]. The 
solution space is not empty since the intertwiner 

s(e u e D )VD~\\N Al N M . . . N Ad > +\N B e BAl - Ao > (3.23) 

is annihilated by M^v, which can be easily checked when using the results of the volume operator 
acting on the solution space of the full set of linear simplicity constraint operators. In order to 
turn the expression into a well defined master constraint operator, we have to square it again 
and to adjust the density weight, leading to 

t 



M N = 4 ^ 5/2 " 2D (Vl^l " N JV Z ) VQ 5/2 ~ 2D N\V T Vi\ ~ NiV 1 ), (3.24) 



which is by construction a self-adjoint operator with non-negative spectrum. We remark that it 
was necessary to use the fourth power of the classical constraint for quantisation, because the 
second power, having the desired property that its solution space is not empty, does not qualify 
as a well defined master constraint operator in the ordering we have chosen. There exists however 
no a priori reason why one should not take into account master constraint operators constructed 
from higher powers of classical constraints [39]. Curiously, the quadratic simplicity constraint 
operators as given above do not annihilate the solution displayed. Clearly, the calculations will 
become much harder as soon as vertices with a valence higher than D are used, since the building 
blocks of the volume operator will not be diagonal on the intertwiner space. 



This type of quantisation is further discussed in section 4.3, where a possible application 
to using EPRL intertwiners is outlined. In contrast to the earlier assumption of D being even 
in order to have an explicit expression nJiyt), we can also perform the mixed quantisation us- 
ing ra 7 nj « (ir aKj l TT a Kj — C 7 } 1 j) f° r Euclidean internal signature £ = 1 and the constraint 
Nj(n n J (it) — N N J ) ~ 0. For the application proposed, we will only need that the correspond- 
ing master constraint can be regularised such that it vanishes when not acting on non-trivial 
vertices, which can be achieved as before. 



4 Comparison to Existing Approaches 

In this section, we are going to comment on the relation of existing results from the spin foam 
literature to the proposals in this paper. In short, the main conclusion will be that in the 
case of four spacetime dimensions, many results from the spin foam literature can be used also 
in the canonical framework. However, they fail to work in higher dimensions due to special 
properties of the four dimensional rotation group which are heavily used in spin foams. We will 
not comment on results based on coherent state techniques [TBI E31 EH [35] since we do not see a 
resemblance to our results which do not make use of coherent states in any way. Nevertheless, 
similarities could be present as the relation between the EPRL |24j and FK [16] models show. 
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4.1 Continuum vs. Discrete Starting Point 



The starting point for introducing the simplicity constraints in the spin foam models is the 
reformulation of general relativity as a BF theory subject to the simplicity constraints, and thus 
similar to the point of view taken in this series of papers. The crucial difference however is 
that while spin foam models start classically from discretised general relativity, the canonical 
approach discussed here starts from its continuum formulation. When looking at the simplicity 
constraints, this difference manifests itself in the choice of {D — l)-surfaces over which the the 
generalised vielbeins (i.e. the bivectors in spin foam models) have to be smeared. Starting from 
a discretisation of spacetime, the set of (D — l)-surfaces is fixed by foliating the discretised space- 
time. Restricting to a simplicial decomposition of a four-dimensional spacetime as an example, 
these would e.g. be the faces / of a tetrahedron t in the boundary of the discretisation. It 
follows that one can take the bivectors B IJ integrated over the individual faces of a tetrahedron, 
By it) '■= ffB IJ , as the basic variables and the quadratic (off)-diagonal simplicity constraints 
read [51] 

C ff := e IJK LB f {t) IJ B f {t) KL diagonal simplicity, V/ G t (4.1) 

C ffl := e IJKL B f (t) IJ B f (t) KL off-diagonal simplicity, V/,/' G t. (4.2) 

In the continuum formulation however, we have to consider all possible (D — l)-surfaces, and 
thus also hypersurfaces containing the vertex v dual to the tetrahedron t. The resulting flux 
operators a priori contain a sum of right invariant vector fields R} e J acting on all the edges e 
connected to v. While this poses no problem for the diagonal simplicity constraints which act 
on edges of the spin networks as shown in [3], the off-diagonal simplicity constraints arising 



when both surfaces contain v are not given by (4.2), but by sums over different Cff, see [3] for 



details. It can however be shown by suitable superpositions of simplicity constraints associated 



to different surfaces that (4.2) is actually implied also by the quadratic simplicity constraints 
arising from a proper regularisation in the canonical framework. This statement is non-trivial 
and had to be proved in [3]. Thus, we can also in the canonical theory consider the individual 



building blocks (2.2) as done in section[2]of this paper. Furthermore, the same is also true when 



using linear simplicity constraints, i.e. the properly regularised linear simplicity constraints in 



the canonical theory imply that all building blocks (3.3) vanish. 

We also note that there is no analogue of the normalisation simplicity constraints [TT] in the 
canonical treatment since the generalised vielbeins do not have timelike tensorial indices after 
being pulled back to the spatial hypersurfaces. 



4.2 Projected Spin Networks 

Projected spin networks were originally introduced in |50| I5T] to describe Lorentz covariant 
formulations of Loop Quantum Gravity, meaning that the internal gauge group is SO (1,3) (or 
SL(2,C)) instead of SU(2). The basic idea is that next to the connection, the time normal 
field, often called i or x in the Spin Foam literature, becomes a multiplication operator since 
it Poisson-commutes classically with the connection. Since the physical degrees of freedom of 
Loop Quantum Gravity formulated in terms of the usual SU(2) connection and its conjugate 
momentum are orthogonal to the time normal field, one performs projections in the spin networks 
from the full gauge group SO(l, 3) to a subgroup stabilising the time normal. Since the projector 
transforms covariantly under SO(l,3), a (gauge invariant) projected spin network is already 
defined by its evaluation for a specific choice of the time normals and the resulting effective 
gauge invariance is only SU(2), which exemplifies the relation to the usual SU(2) formulation in 
the time gauge x 1 '(= N 1 ) = (1,0,0,0). 
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Despite its close relation to the techniques used in this paper and its merits for the four- 
dimensional treatment, there are several problems connected with using this approach in the 
canonical framework discussed in this series of papers which we will explain now. While the 
extension of projected spin networks to different gauge groups has already been discussed in 
|51j . there is a subtle problem associated with the part of the connection which is projected out 
by the projections, that could not have been anticipated by looking at Loop Quantum Gravity 
in terms of the Ashtekar-Barbero variables. There, the physical information in the connection, 
the extrinsic curvature, is located in the rotational components of the connection. To see this, 
consider in four dimensions the 2-parameter family of connections discussed in [2"Q 

A aU = IS + P K *U + l\tU KL KaKL, (4.3) 

where 7 corresponds to the Barbero-Immirzi parameter restricted to four dimensions and (3 is 
the new free parameter appearing in any dimension. K a jj decomposes as [U [2] 

K aIJ = 2N {I K a]J] + K*T + K^T free > (4-4) 

where K a jj means that N 1 K a jj = N J K a jj = and the trace / traceless split is performed 
with respect to the hybrid vielbein. The extrinsic curvature which we need to recover from A a u 
is located in K a j, whereas K^fj e vanishes by the Gaufi constraint and K^fj e free is pure gauge 
from the simplicity gauge transformations. 

Now setting f3 = and N 1 = (1, 0, 0, 0) in four dimensions, we recover the Ashtekar-Barbero 
connection and see that the physical information is located in the rotational components of A a jj. 
It thus makes sense to project onto this subspace in the projected spin network construction, i.e. 
we are not loosing physical information. On the other hand, setting 7 = in four dimensions 
or going to higher dimensions, we see that a projection onto the subspace orthogonal to N 
annihilates the physical components of the connection. This would not be necessarily an issue if 
one would just project the projected spin network at the intertwiners, but when one tries to go to 
fully projected spin networks as proposed in [50]. Then, since one would take a limit of inserting 
projectors at every point of the spin network, the physical information in the connection would 
be completely lost. 

Next to this problem, there are other problems associated to taking an infinite refinement 
limit for projected spin networks as discussed by Alexandrov [50] and Livine |51j . e.g. that fully 
projected spin networks are not spin networks any more (since they only contain vertices and no 
edges) and, connected with this problem, that the trivial bivalent vertex, the Kronecker delta, 
is not an allowed intertwiner. Similar problems have been encountered in section [3j i.e. while 
the vertex simplicity constraints could be solved by a construction very similar to projected spin 
networks where one projects the incoming and outgoing edges at the intertwiner in the direction 
of the time normal N 1 , imposing the linear simplicity constraint on the edges, one would have to 
insert "trivial" bivalent vertices of the form N 1 N J at every point of the spin network, whereas 
one would need to insert the the Kronecker delta 5 IJ to achieve cylindrical consistency while 
maintaining a spin network containing edges and not only vertices. 

Thus, the main problem with using (fully) projected spin networks is connected to the fact 
that we do not know of an analogue of the Barbero-Immirzi parameter in higher dimensions 
which would allow us to put the extrinsic curvature also in the rotational components of the 
connection. In four dimensions on the other hand, this problem would be absent and one would 
be left with the issue of refining the projected spin networks, which is however also present 

7 Note that the definitions of the parameters are different in [2] for calculational simplicity, but here we prefer 
this parametrisation to make our point clear. 
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in section [3] of this paper. Therefore, using projected spin networks in four dimensions with 
non-vanishing Barbero-Immirzi parameter is an option for the canonical framework developed 
in this series of papers and the known issues discussed above should be addressed in further 
research. 



4.3 EPRL Model 

The basic idea of the EPRL model is to implement the diagonal simplicity constraints as usual, 
but to replace the off-diagonal simplicity constraints by linear simplicity constraints which 
are implemented with a master constraint construction [21] or weakly |52| . Furthermore, the 
Barbero-Immirzi parameter is a necessary ingredient. We restrict here to the Euclidean model 
since its group theory is much closer to the connection formulation with compact gauge group 
SO(D + 1). While the diagonal simplicity constraints give the well known relation 

(j + ? = (^)\r) 2 , (4.5) 



.7-1 

the master constraint for the linear constraints gives [23], up to h correction; 



where k is the quantum number associated to the Casimir operator of the SU(2) subgroup 
stabilising N 1 . Depending on the value of the Barbero-Immirzi parameter, either k = j + + j~ 
or k = \ j + — j~\ is selected by this constraint. The EPRL intertwiner for SO(4) spin networks 
with arbitrary valency |29| is then constructed by first coupling the two SU(2) subgroups of 
SO(4) holonomies in the representations calculated along incoming and outgoing edges 

to the intertwiner, to the k representation. Then, the k representations associated to each edge 
are coupled via an SU(2) intertwiner and the complete construction is then projected into the 
set of SO(4) intertwiners. 

An alternative derivation proposed by Ding and Rovelli [52J makes use of weakly implement- 
ing the linear simplicity constraints, i.e. restricting to a subspace T~L ext such that 



C 



V =0 V |0>,|V>e?r xt . (4.7) 



In this approach, one can also show that the volume operator restricted to T~L ext has the same 
spectrum as in the canonical theory, which is an important test to establish a relation between 
the canonical theory and the EPRL model. 

Closely related to what we already observed in the previous subsection on projected spin 
networks, the EPRL model makes heavy use of the fact that SO(4) splits into two SU(2) sub- 
groups and that the Barbero-Immirzi parameter is available in four dimensions. Thus, we would 
have to restrict to four dimensions with non- vanishing 7 if we would want to use EPRL solution 
to the simplicity constraints. One upside of this solution when comparing to our proposition for 
solving the quadratic constraints is that no choice problem occurs, i.e. if we map the quantum 
numbers of the EPRL intertwiners to SU(2) spin networks, a change of recoupling basis in the 
SU(2) spin networks results again in EPRL intertwiners solving the same simplicity constraints. 



8 Note that these ft corrections are necessary since the master constraint, by construction, has the same solution 
space as the original constraint |49| . i.e. C^C\ip) — implies C \ip) = 0. In the master constraint language, one 
subtracts an operator from the master constraint which vanishes in the classical limit to obtain a sufficiently large 
solution space. 
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The problem of stability of the solution space % of the simplicity constraint under the action 
of the Hamiltonian constraint is however, to the best of our knowledge, not circumvented when 
using EPRL intertwiners. 

Also, in order to use the EPRL solution in the canonical framework, one would have to discuss 
exactly what it means to use linear and quadratic simplicity constraints in the same formulation, 
i.e. if one can freely interchange them and how continuity of the time normal field is guaranteed 
at the classical level if one changes from the quadratic constraints to linear constraints from 
one point on the spatial hypersurface to another. The mixed quantisation proposed in section 



3.5 can be seen as an attempt to using both the time normal as an independent variable as 



well as quadratic simplicity constraints. In this case, the main difference is the presence of an 
additional constraint relating the time normal constructed from the generalised vielbeins to the 
independent time normal (which could be used in the linear simplicity constraints). In section 



3.5, this additional constraint was regularised as a master constraint which acts only on vertices. 



Taking the point of view that one can freely change between using the quadratic constraints 
plus this additional constraint or the linear constraints, one could choose the linear constraints 
for vertices and the quadratic constraints for edges. Since we can use a factor ordering for the 
master constraint where a commutator between a holonomy and a volume operator is ordered 
to the right, the master constraint would vanish on edges and only the quadratic simplicity 
constraints would have to be implemented, which are however not problematic. At vertices, we 
would be left with the linear constraints and could use the EPRL intertwiners. Thus, the EPRL 
solution seems to be a viable option in four dimensions. Whether one considers it natural or 
not to use both linear and quadratic constraints in the same formulation is a matter of personal 
taste. Nevertheless, it would be desirable to have only one kind of simplicity constraints. 

A further comment is due on the starting point of spin foam models, which is a BF-theory 
subject to simplicity constraints. It has been argued by Alexandrov [53] that the secondary con- 
straints resulting from the canonical analysis, i.e. the /^-constraints on the torsion of A a j j from 
our companion paper [2j, should be taken into account also in spin foam models. In the present 
canonical formulation, these constraints were removed by the gauge unfixing procedure [2] and 
thus do not have to be taken into account here. The requirement for the validity of this step was 
to modify the Hamiltonian constraint by an additional term quadratic in the /^-constraints 
(the gauge unfixing term) which renders the simplicity constraints stable. While this ensures 
that we have to deal only with the non-commutativity of the (singularly smeared, or quantum) 
simplicity constraints in the present paper, the converse does not necessarily follow: Since the 
Hamiltonian constraint one obtains from the canonical analysis of BF-theory subject to simplic- 
ity constraints, the classical starting point of spin foam models, is not the modified Hamiltonian 
constraint considered here, but the one which results in the secondary ^^-constraints, it does 
not follow that these secondary constraints do not have to be taken into account in spin foam 
models. On the other hand, the present formulation hints that it might be possible to construct 
a spin foam model subject to simplicity constraints (and not D^-constraints) which coincides 
with the dynamics defined by the modified Hamiltonian constraint. In fact, it was recently 
shown that the transfer operator of spin foam models can be written as T =^Wfj, (here for the 
EPRL model) [54] , where fi projects onto the solution space of the simplicity constraint. Taking 
into account the philosophy of spin foam models to impose the simplicity constraint at every 
time step in order to ensure that the second class ^^-constraints are satisfied, it is conceivable 
that the gauge unfixing term of the Hamiltonian constraint in [TJ [2] could emerge from these 
/i-projections when taking the continuum limit of the spin foam transfer operator. Thus, in the 
light of plausible arguments for both sides, only an explicit calculation will be able to decide 
this issue. 

As a last remark, we point out that the non-commutativity of the linear simplicity constraints 
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in the EPRL model results from using 7 / and thus we are not faced with this problem in 
higher dimensions. Essentially, as discussed in more detailed in remark[5]of section 3.4, while the 
rotations stabilising N 1 form an SO(D) subgroup of SO(D + 1), the linear simplicity constraints 
in four dimensions with 7 7^ and (3 7^ do not generate such a subgroup. 



5 Discussion and Conclusions 

Let us briefly discuss the results of this paper and judge the different approaches. 

First, the mechanism for avoiding the non-commutativity in the quadratic simplicity con- 
straints discussed in section[2]is new to the best of our knowledge and we do not see any indication 
that the solution space is identical to previous results (up to the fact that it has the same "size" 
as SU(2) spin networks). In the spin foam literature, the linear simplicity constraints are cor- 
nerstones of the new spin foam models and have been introduced since the quadratic simplicity 
constraints acting on vertices do not commute. While the methods for treating supergravity 
discussed in [5] necessarily need an independent time normal and thus suggest using linear sim- 
plicity constraints, there is no need for the linear constraints in pure gravity (except for the fact 
that they exclude the topological sector in four dimensions). Therefore, one should not dismiss 
the quadratic constraints, especially since the linear constraints come with their own problems 
in the canonical approach. The solution presented in section [2] is certainly not free of problems, 
most prominently the choice of the maximal commuting subset, but its close relation the SU(2) 
based theory and the (natural) unitarity of the intertwiner map to SU(2) intertwiners make it 
look very promising. 

The linear simplicity constraints come with their own set of problems, many of which were 
already known in the spin foam literature. While the results of section [2] would naturally 
lead us to consider the quadratic constraints, the connection formulation of higher dimensional 
supergravity developed in [5] makes it necessary to use an independent time normal as an 
additional phase space variable. This time normal would naturally point towards using linear 
simplicity constraints, although the mixed quantisation of section |3.5| could avoid this. Since 
there is no anomaly appearing when using the linear simplicity constraints (with 7 = in four 
dimensions), we should implement them strongly. However, this leads to a solution space very 
different from the SU(2) spin networks. At this point, it seems to be best to let oneself be 
guided by physical intuition and the results from the quadratic simplicity constraints as well as 
the desired resemblance to SU(2) spin networks. Ad hoc methods for getting close to this goal 
have been discussed in section [3T4| We however stress that these methods are, as said, ad hoc 
and they don't follow from standard quantisation procedures. The mixed quantisation discussed 
at the end of section [3] also does not seem completely satisfactory, especially since the master 
constraint ensuring the equality of the independent normal N 1 and the derived normal n 1 (it) is 



very complicated to solve. Nevertheless, in section 4.3, an application to EPRL intertwiners is 
outlined which could avoid this problem by using linear simplicity constraints for the vertices. 
The strength of the mixed quantisation is thus that it provides a mechanism to incorporate both 
the quadratic simplicity constraints as well as an independent time normal in the same canonical 
framework, which is what is done on the path integral side in the EPRL model. 

A comparison to results from the spin foam literature, especially projected spin networks and 
the EPRL model, shows that many of the problems connected with using the linear simplicity 
constraints have already been known, partly in different guises. While using these known results 
in our framework seems to be a viable option in four dimensions, we are unaware of possible ways 
to extend them also to higher dimensions since main ingredients are a non-vanishing Barbero- 
Immirzi parameter as well as special properties of SO(4). 

In conclusion, we reported on several new ideas of how to treat the simplicity constraints 
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which appear in our connection formulation of general relativity in any dimension D > 3 
[H [21 O HI [5j [6] and found that none of the presented ideas are entirely satisfactory at this 
point and further research on the open questions needs to be conducted. We hope that the 
discussion presented in this paper will be useful for an eventually consistent formulation. 
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A Simple Irreps of SO(D + 1) and Square Integrable Functions 
on the Sphere S D 



There is a natural action of SO(D + 1) on F G % := L2(S D ,d^) given by ir(g)F(N) : = 
F(g~ 1 N). The ir{g) are called quasi-regular representations of SO(D + 1). The generators 
in this representation are of the form tjj = ^(-^jNj — qj^jNi) and are known to satisfy the 
quadratic simplicity constraint t^jt kl -\ = [TT]. These representations are reducible. The 
representation space can be decomposed into spaces of harmonic homogeneous polynomials 
$j D + 1 J of degree I in D + 1 variables, L2(S D ) = Yl'iZo Sj D+1 ' 1 . The restriction of ir{g) to these 
subspaces gives irreducible representations of SO(D + 1) with highest weight A = (Z, 0, 0), 
I G N. These are (up to equivalence) the only irreducible representations of SO(-D + l) satisfying 
the quadratic simplicity constraint [TT] and therefore are mostly called simple representations 
in the Spin Foam community, which we will adopt in this work. Note that these representations 
have been studied quite extensively in the mathematical literature, where they are called most 
degenerate representations [55j [561 E3 > (completely) symmetric representations [561 ESI EHl [60] 
or representations of class one (with respect to a SO(D) subgroup) [33]. The latter is due to the 
fact that these representations of SO(D+l) are the only ones which have in their representations 
space a non-zero vector invariant under a SO(-D) subgroup, which is exactly the definition of 
being of class one w.r.t. a subgroup given in [43] . An orthonormal basis in S) D+1 ' 1 is given by 
generalisations of spherical harmonics to higher dimensions [4.3] which we denote r,f(N), 

f Z?(N) E$(N) dN = 5\,5% (A.l) 

where K denotes an integer sequence K := (Ki, . . . , Kd-2, ±-^d-i) satisfying I > K i > . . . > 
Kd-i > and analogously defined M. Fi{N) £ S) D+l ' 1 can be decomposed as Fi{N) = 
a j?^f (-^0 wnere the sum runs over those integer sequences K allowed by the above in- 
equality. Since L2(S D ) = YlhLo $j D+1,1 , any square integrable function F(N) on the sphere can 
be expanded in a mean-convergent series of the form [43] 

oo 

F ( N ) = Y,Y, al Kn l{N) - (A - 2) 

1=0 jt t 

Consider a recoupling basis [61] for the ONB of the tensor product of irreps: Choose a 
labelling of the irreps Aj, ...,Ajv- Then, consider the ONB 

Ai, A N ; A12, Am, Ai...jv-i; A, Mj , (A.3) 

(with certain restrictions on the values of the intermediate and final highest weights). A basis 
in the intertwiner space is given by 

Ai,...,A^;Ai2,Ai 23 ,...,Ai...jv_i;0,0^ , (A.4) 

(with certain restrictions). A change of recoupling scheme corresponds to a change of basis in 
the intertwiner space. A basis in the intertwiner space of N simple irreps is given by 

Ai,...,A iV ;Ai2 I Ai23,-,Ai... i v-i;0,o\, (A.5) 

(with certain restrictions), since in the tensor product of two simple irreps, non-simple irreps 
appear in general [601 [59] . 

A 2 A 2 -fc 

(A 1 ,0,...,0)®(A 2 ,0,...,0) = 5^ ^(Ai + A 2 -2fc-Z,Z,0,...,0) (A 2 < Ai). (A.6) 

fc=0 /=o 
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